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A LIFT-CANCELLATION TECHNIQUE IN LINEARIZED 
SUPERSONIC- WING THEORY 1 

By Harold Mieels 


SUMMARY 

A lift-cancellation technique is 'presented for determining 
load distributions on thin wings at supersonic speeds. The 
loading on a wing having a prescribed plan form is expressed as 
the loading of a known related vying ( such as a two-dimensional 
or a triangular wing) minus the loading of an appropriate can- 
cellation wing. 

A general expression is derived for the load distribution over a 
cancellation wing. The expression is valid when the plan-form, 
edge (on the cancellation wing) separating a region of zero 
upwash from a region of known loading is everywhere subsoni- 
cally inclined to the free stream. The boundary conditions can 
be satisfied for both subsonic leading and subsonic trailing 
plan-form edges on the prescribed wing. 

The lift-cancellation technique can be used to find the loading 
on a large variety of wings. Applications to swept wings having 
curvilinear plan forms and to wings having reentrant side edges 
are indicated. 

INTRODUCTION 

The method of lift cancellation for obtaining the lift dis- 
tribution on thin wings at supersonic speeds was first sug- 
gested in reference 1. The lift distribution on a given wing 
is determined by canceling excess lift-, through the use of a 
“cancellation wing,” on a related plan form having a known 
loading. This approach has been applied by several authors 
(for example, references 2 to 4). The expressions provided 
in reference 1 are applicable for wings that can be generated 
by the superposition of conical fields. 

A procedure is presented in reference 5 for determining lift 
on a more general class of plan forms than can be handled by 
conical superposition. The method u tiliz es a surface dis- 
tribution of doublets and an inversion by means of Abel's 
integral equation and is equivalent to a lift cancellation. 

This report, prepared at the NACA Lewis laboratory, re- 
tains certain features of reference 5 (that is, the use of a sur- 
face distribution of doublets and an inversion by means of 
Abel's integral equation), whereas other features are simpli- 
fied and generalized. The simplification consists in elimin- 
ating steps in the procedure for obtaining lift distributions. 
The generalization consists in determining a solution that 
can be made to satisfy the boundary conditions for either a 
subsonic leading edge or a subsonic trailing edge (Kutta con- 
dition). The method of reference 5 yields only the Kutta 


solution. The lift-cancellation technique developed herein 
is illustrated by several examples. 

In a concurrent investigation (reference 6), source dis- 
tributions and integral-equation formulations are applied to 
obtain the loading on a special series of cancellation wings. 
Reference 7 employs some of these cancellation wings for the 
determination of lift and moments on swept wings. 

THEORY 

The usual assumptions of an inviscid fluid and small per- 
turbations are made. The velocity field consists of the free- 
stream velocity U (taken in the positive x-direction) plus the 
perturbation velocities u, v, and w. The wing boundary 
conditions are specified in the z=0 plane. 

The local lift coefficient AC, may be expressed in terms of 
Au; that is, 

A/~r Vb Vt 2 (u T -u B ) 2Au 

AC P — v jj~ ( 1 ) 

(All symbols used in this report are defined in appendix A.) 
Inasmuch as the local lift coefficient is directly proportional 
to Au, Au will be referred to as “lift” in later developments. 

LIFT-CANCELLATION METHOD 

The lift distribution on a given wing is to be determined by 
canceling excess lift on a related wing with a known loading. 
The method is illustrated in figure 1 . The wing for which the 
lift distribution is desired is shown in figure 1 (a). The solu- 
tion can be expressed as the two-dimensional wing (fig. 1 (b)) 
minus a cancellation wing (fig. 1 (c) ) . The loading in region I 
of the cancellation wing equals the loading in the correspond- 
ing region of the two-dimensional wing and the upwash w 
in region II of the cancellation wing is zero. The two- 
dimensional wing minus the cancellation wing satisfies the 
boundary conditions for the flow about the given wing and 
is the desired solution. 

The fundamental problem in the lift-cancellation method 
is then to determine the lift in region II of a cancellation 
wing subject to the condition w=0 in this region and with 
the assumption of a known loading in region I. Solution of 
this problem is presented in the following sections. 

DERIVATION OF LIFT-CANCELLATION EQUATIONS 

The lift distribution in region II will be expressed in terms 
of quantities in region I. 


•Supersedes NACA TN 214J, “Lift-Cancellation Technique In Linearized Supersonic- Wing Theory'’ by Harold MIrels, 1SKL 


213637 — 53 6 


65 



66 


REPORT 1004 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


Consider the cancellation wing shown in figure 2. The 
portion of the leading edge to the left of the origin coincides 
with a Mach line. The portion of the leading edge to the 
right (designated r=rj(s)) is shown as a supersonic edge, 
although no restrictions as to a subsonic or supersonic edge 
are imposed. (A plan-form edge is subsonic or supersonic 
depending on whether the component of .the free stream 
normal to the edge is subsonic or supersonic.) The line 
designated r=r 3 (s) separates region I and region II and is 


i 





(O) 

(a) Given wing. 

(b) Two-dimensional wing, 

(o) Cancellation wing. 

Fiouhb L— Superposition to obtain lilt on given wing by canceling lilt on two-dimensional 
wing. (Given wing equals two-dimensional wing minus cancellation wing.) 


assumed to be subsonically inclined to the free stream at all 
points. This line corresponds to a plan-form edge of the 
wing for which the lift distribution is desired. 

General solution for Av on cancellation wing. — The up- 
wash field in the 0=0 plane (associated with an arbitrary 
distribution of vortidty A u and Av) may be written, from 
reference 8, 


j3 2 | f f [(y— y»)Ag + (ig — x t )Au]dx f dyl 

w 2w| JJ T [(x-x o r-p\y-y 0 )V a {) 

The symbol | designates the finite part of an infinite inte- 
gral, as defined in reference 9. Application of the finite-part 


concept to linearized supersonic-wing theory and the evalua- 
tion of the finite part of an infinite integral are discussed in 
references 8 and 10. For the present, it will suffice to state 
the fundamental definition of the finite part of an integral 
with a 3/2-power singularity, namely, 


| f x jbo)dxa f x tf(x 0 )-f(x)}dx 0 2/(a:) m 

LUx-^-J. (x-Xo) 3li ~ ( s-a) 1 '* W 

By a transformation to the Mach coordinates of reference 1 1, 


(8+r) 


V=M (s_r) 


s=||(z+0y) 

28 

elemental area=-^ dr da 


(4) 


equation (2) becomes 
M 


v>— 


8 : 


rr _ (s - s#) ^: +(r ro) brj dr ° ds ° (5) 
J I ' \(r~r„) ( )F 


Upon substitution of the limits of integration, as indicated 
in figure 2, 


w—- 


M< 

8w< 


dA <f 


dr 0 


r ds 0 f Ar, 

)o (s-s^Vnc*.) {r-r 0 )d* 


r* ds,, 

r 

Jo (s-s 0 ) ]/J j 

nc».) (r—r 0 ) z/1 ) 


( 6 ) 


U 



Integrating by parts, noting that A<p=0 at r e ,=ri(s 0 ), and 
recalling the definition of the finite part (equation (3)) yield 


dA ip 


dr 0 


C r ° =Uml_ ^L-S\ r ° C r " A<pdr Q ) 

Jr lW (r— r<,) 1/s -nUr-t^JMO 2 Jr l(0 (r-r,H 

( 7 ) 


s _i r 

2 Jr,(«,) (r— 


A <pdr 0 


r,) 3 ' 2 
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Thus, 


5Ay , 

r« di , r *>r 0 ar ' _ ll r d&g r a ydr„ 

0 (a— $») 3/ * Jr,(0 (r—r 0 ) m 2[ J 0 (*— O'— r tf ) s/s 

( 8 ) 


Similarly, reversing the order of integration (with appro- 
priate changes in limits of integration) , integrating by parts, 
and then returning to the original order of integration 
establish the identity 


dA £ , 

r* ds t c r •_ iff* 5 ^ r a ^*7 

Jo (s-s .) I/2 Jrtft.) (r— 2 |Jo (s-s„) 3/i Jr lW (r-r a ) 3 ' 2 

( 9 ) 


The right side of equations (8) and (9) are identical. Equa- 
tion (6) can now be written as 


M 

f* ds 0 | 

C r A <pdr a 

8x , 

Jo (s-So) 3II J 

n (>,) (r—r,) zn 


( 10 ) 


For points in region II, w = 0 and equation (10) becomes 


or 


where 


r* ds 0 

C r A <pdr a 

Jo (s-s a ) 3 '\ 

Jr itto (r— r „) 3/2 


f G(r,s 0 )ds 0 

Jo (s-s») 3/s 


G(r, s,) 


= C r Ay 
JriC».) O’— 


A <pdr 0 


r„) ,/2 


(Ha) 

(Hb) 

( 12 ) 


Equation (lib) is an integral equation for the unknown func- 
tion G(r,8,). The solution (appendix B) is 


Thus, 


G(r,s 0 )= 0 


(13) 




(a) Region I Intersected by right forward Much line from (r,a). 

(b) Region I intersected by left forward Mach line from (rj) . 


Fiona* 3. — Geometric interpretation of terms in equations (17a) and (17b). 


General solution for load distribution on cancellation 
wing. — The load distribution in region II can be expressed as 

dA <p a br dA ip a da 

Uu ~ fa ~~ br b x bs bx 


= m /&_ _e>\ 

2/S 


A(f>n 


or 


r A <pdr 0 _ Ayidr„ f r Aip u dr 0 

Jn(*.) (r — r«) 3/ * Jr l( ,») (r-r 0 ) 3/J+ Jr j( ,.) (r— r,) 5/ * 


rp A <p u dr a Ay r dr 0 

IJra<».) O'— r a ) t/2 JnCr.) (r— r ,} 3/2 


(14) 


The right side of equation (14) will be considered known. 
Equation (14) is then an integral equation for A <pu. The 
solution (appendix B) is 


or, from equation (15), 



(b b\V Vr— r t (s ) A<pidr e 1 

\br + WL x 'Jr lW (r-r.)Vr a (*)-rJ (16) 


Differentiation yields (see appendix C) 

. . ->Jr~ r t (s) f r * w A«i dr 0 

*■ (r— r 0 )-v/r 2 (s)— r. 


\,„ Vr— r a (s) f r » w Atpidr, 

r JnW (r— r<,)Vr 2 (s)— r 0 


(15) 


Equation (15) indicates that the doublet strength in region 
II, namely Apn, can be obtained by a line integration 
along s„=s in region I. The geometric interpretation of the 
various terms in equation (15) is shown in figure 3 (a). 

It can be shown, by expanding Ai pi about r 0 =T t (s), that 
equation (15) yields a continuous solution (A£n=Api) 
at T=r t (s). (A discontinuity in A <p implies a lifting line 
(reference 12) and is unrealistic.) 


, dr 2 (s) 

da (jSAtfj— ApQdr. 

2/5xVr— r 2 (s) Jn« V r 2 (s)— 

Equation (17a) is the desired expression for the lift distribu- 
tion in region II in terms of quantities in region I. 

Consider A u n to consist of two components, Au u ' and 
Attn", where Att n ' and Attn" are the first and second terms 
on the right side of equation (17a), respectively. Investiga- 
tion of the integrals indicates that at r=r t (s ), Au n '=Aui; 
whereas Au u ", in general, has a half-order singularity. 
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When region II is to the right of region I (fig. 3 (b)), the 
integration for A u n is conducted along the line r a =r and 
may be written as 

A ^ f* a(r) kuidsp 
U v J* lW (« — So) -\ 8 i(r) — So 

1 dSjjr) 

dr paM (j3Aui-|-A«i)ds a ^ 

2j3 ttV® — S j(r) *J$i(r)—s c ^ ^ 

Discussion of equations (17a) and (17b). — In the par- 
agraph preceding equation (2), the line r=r 2 (s) was described 
as subsonically inclined at all points to the free stream. 
This condition is necessary so that the inner integral in equa- 
tion (11a) (that is, 6(r,s 0 )) can be equated to zero for all 
points in region II. If this restriction on r=r 2 (s) is not 
satisfied, the development beyond equation (11a) becomes 



& r x 

A«i =■=— Auidx, (18a) 

The integration is conducted along lines of constant y. 
In the second case (fig. 4(b)), region II is upstream of region I 
(along the line r=r 2 (s)) and the expression for Ai>i in region I 
(for y< 0) is 

or A “ nii ' + £,> <i8b) 

Equation (18b) indicates that a knowledge of Au u is required 
in order to find Avz. But Ai^ must be known (equa- 
tion (17a)) before Au u can be found. Thus, the solution for 
Aun from a specified Au t is not unique for the configuration 
of figure 4(b) and an additional boundary condition must 
be imposed. The line r=r 2 (s), however, corresponds to a 
plan-form edge of the airfoil whose load distribution is 
desired. The situation indicated in figure 4(b) occurs when 
r=r 2 (s) corresponds to a subsonic trailing edge. The 
additional condition to be imposed is therefore the Kutta 
condition. In terms of the cancellation wing, this condition 
requires that the perturbation velocities be continuous in 
crossing r=r 2 (s). 

Solution for A u n satisfying Kutta condition at r=r 2 (s). — 
It will now be shown that when the Kutta condition is 
imposed at r=r 2 (s), the appropriate Ai^ distribution is such 
as to make the second integral in equation (17a) identically 
zero; that is, 



(a) Region I upstream of region II (along r-rs(*)). 

(b) Region II upstream of region I (along r-ra{*)). 


Fiouee 4.— Possible relations between regions I and II In regard to determination of dri. 


r «> m =^ ir , =0 

Jn(.) [r 2 (s)— r 0 ] l/2 


or, inasmuch as fiAui—Av-i—M - A<pl ; 

Or o 


j: 


5Ag>i 
r » w dr 0 


dr 0 


(«) [r 2 (s)— r,] l/2 


This concept and its proof foDow from a suggestion of 
H. S. Ribner of the NACA Lewis laboratory. 

Thus, from equations (7), (12), and (13), 


X 


dA (p 
dr. 


dr 0 


r,t.) (r-r,) 1 '* 


(19) 


invalid. The derivation of cancellation equations when 
r—r 2 (s) is supersonic was not undertaken because such 
problems can be solved more simply by other methods. 

In regard to the boundary conditions, it has been assumed 
that AUx is specified. Equations (17a) and (17b), however, 
indicate that a knowledge of A?^ is also required in order 
to obtain a solution for Au n - With regard to the deter- 
mination of Avi, two possibilities exist, as illustrated in figure 4. 
In the first case (fig, 4(a)), region I is upstream of region II 
(along the line r=r 2 (s)) and A^ is uniquely defined by the 
specified Aiij according to the relation 


for all points (r,s) in region II. Therefore, 


dAgu , vntpn , 

&r. dr ° C r ?>r, dr ° 
Jr lW (r-r 0 ) in ~ J r a w (r-r,) 1 ' 2 


dA^n 


When the limit as r approaches r 2 (s) is taken, equation (20) 
becomes 


bA<p! J* 
dr, d ° 

— li tn 

dr ~ 

C r dr, aTo 

[J , 3 (s)-r 0 ] 1 '' 2 

r ->r 2 («) 

L Jr 2 C»J (r— r„) 1/3 J 


( 21 ) 
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However, 3A tpifbr, must be continuous in tbe vicinity of 
r 2 (s). (The perturbation velocities on the basic wing can be 
discontinuous only along Mach lines or along plan-form edges. 
Inasmuch as r=r 2 (s) is neither of these cases, all derivatives 
of A^r must be continuous in the vicinity of r=r 2 (s).) 
When the Kutta condition is imposed, dA^n/dr, is there- 
fore also continuous (and bounded) in the neighborhood of 
r—r s (s). Then with the use of a mean value for bA<p n fdr„, 


Tim 

MliW 


5A ipu. 


= lim 


i r 

L Jr 2 (*) (r— r a ) I/2 J 

r/bA^uN r _dr. 1 =0 

L\ Sr 0 / (r,(i)o.<7) Jr,« (r— r.) 1/2 J 


Therefore, 


X 


bA<px 


dr. 


ri(i) [r 2 (s)— r .] 1 ' 2 


(22) 


(23) 


which tfas to be proved. 

The solution for Au n that satisfies the Kutta* - condition at 
r=rJ(s) is then, from equations (17a) and (23), 


Ait 


n- 


Vr— r 2 (s) f r 2 (,) hu r dr, 


JnW (r— ; 


(24a) 


-r 0 )^r 2 (8)—r 0 

for the wing of figure 3 (a) . Similarly, 

. Ys-s,M Auids a 

A u u =- 7 . I , . (24b) 

ir JnW (s— s 0 )Vs2(r)— s„ 

for the wing of figure 3 (b) . 

An alternative derivation of equations (24a) and (24b) 
(appendix D) indicates that only solutions satisfying the 
Kutta condition will result from the integral equation for- 
mulation of reference 5. 

Sidewash in region n. — An expression for Acn can be ob- 
tained by differentiating equation (15) with respect to y. 
The result is 


Ari 


v'r— r 2 («) Avi dr 


rr 2 W 

JnC») (r — j 


dr 2 (s) 


(r— r 0 )Vr 2 (s)— r ( 
ds C r * l,) igAu r — Aa r 

2 tt, r— r 2 (s) JnM Yris (s)— r 0 
Similarly, for region II to the right of region I, 


dr 0 


and 


i/r — r 2 (s) r r * W Ari dr „ 

Vu ~ ir Jnw (r — r 0 )-\Jr 2 (s) — r„ 

-yjs— s 2 (r) Aigcfsp 

rrc ~ W J«:W (s— S 9 )Vs 2 (r)— s. 


It should be noted that when r=r 2 (s) corresponds to a 
subsonic trailing edge, Ar I; as well as Arn, is not generally 
known. The preceding expressions are therefore p rimar ily 
useful for those problems where r—r 2 (s) corresponds to a 
subsonic leading edge. 


APPLICATIONS 

The loading in region II of a cancellation wing is given 
by the line integrals of equation (17a) or (17b). When the , 
Kutta condition is imposed at a subsonic trailing edge, the 
expressions reduce to equations (24a) and (24b). These 
equations can be used to find the load distribution on a large 
variety of wings. Wings with curvilinear plan forms or 
arbitrary camber are examples. In each case, however, 
the solution for the related wing must be known. 

The equations are applied in several illustrative examples. 
Only the solution associated with the cancellation wing is 
considered. The complete solution consists of the loading 
of the related wing minus the loading of the cancellation wing. 

LEADING-EDGE AND SIDE-EDGE CANCELLATIONS 

In leading-edge and side-edge cancellations, the lift to be 
canceled is upstream or to the side of the plan form for which 
the loading is desired (figs. 5 and 6 ). 

Tip region of swept wing. — The loading in the region 
influenced by the side edge (H„ and IK of fig. 5) of a swept 
wing having a subsonic leading and a supersonic tra iling edge 
can be obtained by canceling excess lift on a triangular wing. 
The Kutta condition is applied across the portion of r=r 2 (s) 
influencing region II 6 . The lift to be canceled in region I is 
(reference 13, equation (23)) 

Am— r gfc — - r =gg£±a (25) 

Ye 2 * 2 - ,8V Y®*(® +r) 2 — (a— r)* 

where H and 8 are constants defined in appendix A. The 
doublet distribution in region I„, again from reference 13, is 

A< Pl =H-yjei t x' i — Y^s+r) 2 — (s— r) 1 


Apii= Ys = £2W P 2(rt Av^ds, t 

T " *l(r) (s—8 0 )^S i (r)—S e 
, , ds 2 (,r) 

_ ds gAgi+Agi j. 

2xYs— s 2 (r) Y«i (»•)—«. 

When the Kutta condition applies, these equations become, 
respectively, 


from which 


a,, _ £)Ay>I - _ —Ep*y —Hpjs—r) 

*■ Yfl***-flV Y^(«+r) 2 -(s-r) s 


(26) 


The sidewash distribution in region I 6 (that is, Av^ could 
be found by an integration of the type indicated in equa- 
tion (18b) . A knowledge of Avi t , however, is unnecessary in 
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the present problem because the Kutta condition is applied for 
region II>. 

The loading in region II a is obtained by substituting 
equations (25) and (26) , with r replaced by r 0 , into equation 
(17 a), which yields 


Attn 


Hd 2 -\lr—r 2 (s) f r aW 


H\ 1 


V 

Hf], 

2 ttV? 1 — r 3 (s) JnW V r 2 (s) — r„ V 0 s (s + r 0 ) 2 — (s — r 0 ) 


paW 

JnW (r— 


(s+r c )dr 0 


V*W 


(r-r^V^s)— r^V^Cs+^-fs— n) s 

[fffc+r^+fc— r$\dr 0 


(27) 


r 



Fiqubi I.— Cancellation Cor obtaining loading in tip region of swept wing having supersonic 

trailing edges. 

For region II 6 , the Kutta condition applies and 

(s+r<,)dr 0 


Aw : 


g0 a Vr— r 2 (s) 

* JnW (r— 


(r—r,)V> , a(s)—nV^(«+^) ! — («—'■»)* 

(28) 


Equations (27) and (28) reduce to elliptic integrals of the 
first, second, and third kinds upon transforming the variable 
of integration from r„ to w„ according to the relation 


[(1 — 0)8 + a3«o*] 


where 


4 1 + 0 

o j = (1 + 0)r a (s) —(l — 0)g 
Equations (27) and (28) may then be written 


(29) 





HK?*)- 

(30) 


and 


AUn t > 


g0»Vr-r,(«) ((s+r) (1+0) 


TT-jsd 


| (8+rKl+0) ^ nj £)j _ 


[F(|,fc)-F(0,fc)]J 


where 



a=(l+0)r— (l — 0)s 

* =sin "VI 

a 1 =(l + 0)r,(s)— (1 — 0)8 

n— — — 
a 

a 2 =(l+0)r 2 (s)— (1 — 0)s 


Reentrant side edge. — A plan form has a reentrant edge 
if a line of constant y intersects the plan form at more than 
two points. 

The load distribution in the region influenced by the re- 
entrant side edge is to be determined for the wing (unshaded 
region) of figure 6. The side edge is first, for simplicity, the 
straight line r=K 2 s, which is a subsonic trailing edge 
across which the Kutta condition is applied. The side edge 
then alternately becomes a subsonic leading and a subsonic 
trailing edge. The load distribution in region I is simply 

the Ackeret value Au I =^^~, and Av In =0. Regions II*, II*, 

k p 

and !!<■ are considered separately. 



Fioubi 6. — Cancellation for obtaining loading in region influenced by reentrant side edge* 

Region II„: 

2aU 

From equation (24a) with A«i =— > 


AUr 


■\jr—KjS r K *‘ 2 all dr, 


r*2> 2 

J-* B(r— 


d(r r ff ) ■yjKiS T 0 


+ 1 )* 


(32a) 


4 a u ± mr- 

=~z — tan \ - T T 

jSt V r—KtS 

or, in x,y coordinates 

‘“W g (£«-») (32b) 

Region H*: 

A knowledge of Avj 6 is required. From equation (18b) , 

dR aU C ltm * . I z°+Py . 2 aU'r* 

tan ' V fKmtXo-y)**'* p JtftH**} 


2aU 

Vf 


=2 a U. 


lKt+l 


(33) 


(31) 
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The load distribution in region Us is then, recalling that Ar^=0, 


A-Utt = 


__ y r—r 2 (s) 2aUdr, 

r J-* P(r— r a )^r a (s)— r B 


dr 2 (s) ~ j 

J f f* 2aUdr„ f r * w / 2aU—2*U^lE*iE£) . dr ° .1 

2/3 x\/r— r 2 (s) [_J -* -Jr 2 (s)—r„ J* \ \ K 2 — 1/ V r s(«) — ^sj 


Region II C : 

Because the Kutta condition is applied, 


. 4oE7 . . 

Aun o =-^- tan 


_i / r 2 (a)+s 

V r— r s (s) 


TRAILING -EDGE CANCELLATION 


The calculation of lift distributions on swept wings having 
subsonic trailing edges requires cancellation wings of the 
type shown in figure 7. These wings cancel that part of the 
lift of the basic tri ang ular wing that is downstream of the 
trailing edge of the swept wing (references 3, 4, 6, and 7). 
The lift is specified in region I. The lift in regions II and HI 
is to be determined subject to the conditions that w=Q and 
that the Kutta condition applies at r=r t (a) and r=r 2 (s). 


The wing of figure 7 (a) differs from the previously dis- 
cussed cases in that two unknown regions (II and III) are_ 
continuously interacting. A special treatment is required 
in order to obtain the loading in regions H and HL (See, 
for example, reference 6.) Approximate solutions can be 
obtained, however, using equations (24a) and (24b). For 
example, if the load at (r,s) in region II is desired, first 
assume that A« m is known. Then, 


Vr— r 2 (s) r p lW) Aumdr 0 f r * c *> A Ujdr, 1 

la ~ v (r— r a )Vr 2 (s)— r 0 JnW (r—r a )-y/r 2 (s)—r a 

_ 20 J 


An expression for Au ra is, by integration along lines of constant r„ (fig. 7 (a)) , 


Au rn = ^' S Sl ^ r °) 


An approximate expression for Au m is then 


“ r»2(n) 

J ^<v (s— 
. 20 


Au n ds„ ^ p t Cr„) Au z ds 0 

S.)V«l(r.)— S„ JiaCr.) (s— S B )V»lW- 


V’s — Si(r„) Au r ds 0 

V J>2 ('.) (s — s 0 ) V*i(Ja 


v J>2 <x.) (a— s 0 )v s iOv)— s g 

Equation (36), which may now be written in terms of Au r by substituting equation (37) for A« m , becomes 


. rr—r 2 (s) ( 


( p lC,) in r y.s— gj(r 0 ) pifrJ Auids c p» (t) Aui dr 0 ) 

Pi£± (r—r 0 )^r s (s)—r 0 L * ( g — sjy's&o)— sj J r (r— r 0 )Vr a (g)— r B j 

v 2fl 


The first term on the right side of equation (38) approximates 
the contribution of region III to the loading in region II. 
This term, as indicated in reference 6, is negligible for 
the commonly encountered Au z distributions (corresponding 

to steady lift, roll, or pitch) and For those 

cases, equation (38) simplifies to 

a _.\r-r 2 (s) C r * w Auj.drp ,„ Q , 


) p«(«> Auidr s 
Jn(«) (r—r a )^r 2 (s)—i 


The Kutta condition at r=r 2 (s) is satisfied by both equations 


(38) and (39). Equation (39) can be reduced to elliptic inte- 
grals in canonical form by the substitution indicated in 
equation (29). The elliptic integrals may be avoided by 
expanding Au z in a Taylor’s series about r 0 =r 2 {s). Equa- 
tion (39) can then be readily integrated, term by term, to 
yield a very good approximation for A u n in terms of alge- 
braic and trigonometric functions. 

When region I has a partly supersonic leading edge 
(fig. 7 (b)), it is possible to write exact expressions for the 
linearized load distribution in regions IE and HE. For 
example, the load at point (r,s) of figure 7 (b) is 


A tl n = 


\r— r t (s) f f r i« 


U ’nW dr a j~ V-8— Si(r<,) ptfr.) Au z ds a ~ | Aujdr, ) 

e O’ — r„) V r 2 (s) —r 0 L T J*iW (s— s a ) VSi(r,)— s 0 J JnW (r—r 0 )^r 2 (s)—rj 
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which reduces to equation (39) when the contribution of 
region III is neglected. 

From the development of reference 6 it may be concluded 
that equation (39) is sufficiently accurate for most problems 
involving trailing-edge cancellations. The evaluation of 
equation (39) is generally simplified by the expansion of 
Aui in a Taylor’s series about r 0 =r 2 {a). 




(a) Regions II and III continuously Interacting. 

(b) Regions n and m not continuously Interacting. 

FinUBE 7.— Typical cancellation wings lor canceling lift downstream ot subsonic trailing 
edge of swept wings. 

SUCCESSIVE CANCELLATIONS 

A cancellation wing may induce lift that itself must be 
canceled in order to satisfy boundary conditions completely. 
Thus, in figure 8, the cancellation of lift in region I induces 
lift in region I'. The cancellation of lift in region I' induces 
lift in region I", and so forth. Each of these cancellations 


is handled as previously described. These computations are 
very tedious when lift is induced upstream of a subsonic 
leading edge (for example, region I' of fig. 8) , inasmuch as a 
knowledge of the sidewash (Aiv), as well as of the lift dis- 
tribution (Aur), is needed in order to continue the cancella- 
tion process. Numerical methods are generally required. 

Successive cancellations are discussed more extensively in 
references 3 and 4. 



SUMMARY OF ANALYSIS AND APPLICATIONS 

A general expression was determined for the lift distribu- 
tion over a cancellation wing. The expression is valid when 
the plan-form boundary (on cancellation wing) separating 
the region of zero upwash from the region for which the lift 
is specified is everywhere subsonically inclined to the free 
stream. This expression permits the determination of lift 
distributions on a large variety of wings. The boundary 
conditions for either the flow about a subsonic leading edge 
or a subsonic trailing edge can be satisfied. 

The lift-cancellation technique was illustrated for swept 
wings having curvilinear plan forms. Leading-edge, side- 
edge, and trailing-edge cancellations were considered. In 
addition, the loading in a region influenced by a reentrant 
side edge was found. 

Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, January 16, 1950. 
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The 
a = 
a-i 
a 2 

AC, 

c 

Ct 

Efak) 


following symbols are used, in this report: 
= (l+^)r— (1 — fl)s 

= (l + flWs) — (1 — 0)8 

= ( 1 4 - fl)r 2 (s) — (1 —ff)s 

local lift coefficient, — — — 

2 

constant 

root chord of swept wing 

elliptic integral of second kind, 


F{cj>,k) 

C(r,s 0 ) 
H = 

K 

h 

M 

m 

n 

V 


r,r„) 

* ,*•) 


r, 

s, 

U 

u,v,w 


Au 

Av 


APPENDIX A 

SYMBOLS 


P = 

5 

8 = 

M4>,n,k) 


ntegral of 


-\fi — k 2 cc„ 2 j 

u 1 ««» 

VI— Wo 


elliptic integral of first kind, 

*** dug 

So V(1 — 0 > a *)Q.—k*U,*) 

function of r and s t defined by equation (12) 
2aU 


^(ivu^) 


slope of plan-form edge in r,s coordinates, drjds 
modulus of elliptic integrals 
Mach number 

slope of plan-form edge in x,y coordinates, dyfdx 
parameter of elliptic integral of third kind 
local static pressure 

3'^ 

Mach coordinate system (equation (4)) 
free-stream velocity 

perturbation velocities in x-, y~, and z-directions, 
respectively 

u T — u B (proportional to local lift) 
v T —v B 

Cartesian coordinate system 


APPENDIX B 

SOLUTION OF INTEGRAL EQUATIONS 
Consider the following integral equation (in the notation 


angle of attack 

VM 1 — 1 

semivertex angle of tri ang ular wing 
/9 tan 8 

elliptic integral of third kind, 

/"Ski <j> d(j3 0 

Jo 


(1 +716J9 2 ) V(1 —£*«**)( 1 —co, 2 ) 


p 

density 

T 

area of integration 


amplitude of elliptic integrals 

<P 

perturbation velocity potential 

Aip 

doublet strength, pt-vb 

O), 

Regions: 

integration variable 

I 

region on cancellation wing for which loading is 
specified 

Iffjhj ... 

subdivisions of region I 

n 

region on cancellation wing for which w=0 * 

IX,, n*, . * . 

subdivisions of region II 

in 

additional region on cancellation wing for which 
w— 0 

Special designations: 

r=ri(s) 

r as function of s along plan-form boundary 1 

s—8i(r) 

s as function of r along plan-form boundary 1 

y=y i(») 

y as function of x along plan-form boundary 1 

r=z 1 (?/) 

x as function of y along plan-form boundary 1 

r=rj(s) 

r as function of s along plan-form boundary 2 

s=Si(r) 

a as function of r along plan-form boundary 2 

and so forth. 
Subscripts: 

1,2,3 

refers to plan-form boundaries 1, 2, and 3, re- 
spectively 

I, II 

refers to regions I and II, respectively 

B 

bottom surface of z=0 plane 

T 

top surface of z= 0 plane 

0 

variable of integration 


of the appendix in reference 5), where the function f(x) is 
assumed known and the function u(£) is to be dete rmin ed: 


*Hjr^ 


=/. 


[u(Z)—u(x)]dZ 2u{x) 


(Bl) 


(x — £) s/2 (x—a) U2 

After an integration by parts, equation (Bl) may be written 
. /(*) "M (B2) 


2 J* (a;-*) 1 ' 2 

Equation (B2) is now an integral equation of the Abel type. 
The continuous solution for «($) is (reference 14) 


/ \ If* f(x)dx 

u U— 


CB3) 


evaluated at z=f. This result is presented in reference 5. 

Equation (lib) corresponds to equation (Bl) with 
u(Q=G(r, s,) and f(x)=0. The solution for Q(r, «,), accord- 
ing to equation (B3), is then 

Q(r,s 0 )= 0 (13) 

Equation (14) corresponds to equation (Bl) with 

£=r„ a=r 2 (s 0 ) 

: dr~ 


x=r /(*)=_ ** 

Jn(».) (r— i 

u(jO—A<pu 


r,)*' 2 
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The solution for A<p n according to equation (B3) is then 
_ J_ p dr An dr 

2xJ r a (*„) r x (» 4 ) 


> Vz— rJnW ( r “ r - 

Reversing the order of integration and integrating yield 
_ ^z —rjjs,, ) J r M 


(B4) 


A<pi 


A <px dr 0 


M*») (z—r 0 )- s lr 2 (8 l >)—r 0 


(Bfi) 


Equation (B5), evaluated at z=r and s 0 =s, yields 


A<pa z 


Jr 


f r * (,) 

rc J r iM 


A(pi dr 0 


nW {r-r e )-yjr 2 (s)—r 0 


( 15 ) 


The derivation of equation (15) is similiar to that for equa- 
tion (16) of reference 5. 


APPENDIX C 

DIFFERENTIATION TO OBTAIN Au n 


The differentiation indicated in equation (16) 

2/3 /a , b \ f Jr— r *( s ) C ra<,) A<Pidr 0 ~| ,, 
M Un \br ds/j_ r JrjW (r— r 0 )-vVi(s)— rj 

is to be conducted. 

First, 

2)3 


dr 2 (s) 


ds paW Atpzdr, 

M Ul1 2irdr— r 2 (s) Jn<*) (r— r^V^aO 


Vr— r a (s) /b , b\ f 
w ydr^dsjjr 


{r—r 0 )^rt(s)—r 0 
r aW A^>i dr 0 


r iW (r—r 0 )jr s (s)—r 0 
Inasmuch as A^i is a function of r„ and s, 

Ac?! dr 0 _ A^j dr 0 


(Cl) 


Af 

brjr. 


K-X 


ri(«) (r-r,) ! ^!(s)-r, 


and 


(C2a) 


b pa<*> 

(r— 


A<pidr 0 


H 


/ dAy A 
faW \ 3« / 


dr„ 


(r — r„) ^r 2 {s)—r 0 J n W (r — r 0 ) V r a(s)—r 0 

( 1 dr a (s) p« Atpjdr 0 . 

ji“)( 2 ds J ri w (r-r 0 )[r 2 (s)— r„] 87S 


[(AyQr.-riW] dr i(s) 


[(^yi) r »~ r a(»)] dr 2 (8 ) ) 

( 'r—r<>)i]rt(s)—ro ds ) [r-ri(«)]Vr a («)— Fi(«) 


(C2b) 


However, [(A^ I )r.-r 1 c»)] = 0 and, by integration by parts, 


r 2 («) 


A^i dr 0 


! "" 2 X 


bAy 

r a W br, 


1 1 

1 r„)_ 


dr„ 


L (r-r*)[r 2 (s)-r 9 ] 3/2 “Jn(.) (r-r 0 )Vrj(s)-r, 

lim ( 2Ayi T 1 

r*-*ri« ( (r — r„) yr 2 (s) — r 0 J rr (*) ) 
so that equation (C2b) can be written as 


b^ pa(«) Atpidr 0 

dsjr,(») (r—r 0 )^r t {8) 

r&A^i, 

dr 2 (s) p*W L br„ Ar 

ds JrjW (r— r 0 )Vr 2 (s)— r„ 


rr 2 (>) 
—r„ JnW | 


bA <pi 
bs 


- dr„ 


nW (r— r 0 ) Vr a (s)— r 0 


raA^x_ ] __A f p^l d 
-t 8 )J 


(C3) 


Equations (C2a) and (C3) are substituted into equation (Cl) 
and the integrals containing A ^ are integrated by parts, re- 
ducing equation (Cl) to 


Au n = 


■yjr—r 2 (s) pa (,) Auidr„ 


'.I! 


r i<»> (r— r 0 )Vr a (s)— r 0 


dr a (s) 


ds r 
Ir—rds) J r : 


r a(*) (pAuj—Avi)dr 0 


2/3x\V— J*s(s) D r i(«) V r j( s ) —r » 


(17a) 


APPENDIX D 

ALTERNATIVE DERIVATION OF SOLUTION SATISFYING KUTTA CONDITION AT r=r,(*) 


The integral equation formulation in terms of Ap (equa- 
tion (14)) resulted in a solution that was continuous in A <p 
(equation (15)) but discontinuous, in general, in the deriva- 

tive A^=Am (equation (17a)) at r=r a (s). In order to obtain 

a solution continuous in A u, an integral equation may be 
formulated that is similar to equation (14) but in terms of 
Au rather than A<p. The inversion shown in appendix B 
should result in a solution that is continuous in Am but dis- 
continuous in the derivatives of Am at r=r 2 (s). 


Consider equation (10) for the w distribution in the 2=0 
plane. This equation will be differentiated with respect to 
x using a technique introduced in reference 15 (equations (1) 
to (3) therein). The expression for w at any point (r,s) is, 
from equation (10), 


Ad (* {* Atpdrod&o /tiix 
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where r is the area abc in figure 9 (a). The wing is moved 
upstream a distance dx (fig. 9 (b)), keeping the coordinate 
system fixed in space. The expression for the upwash at 
(r,s) now becomes 


. bw , M 


1 1 r r ( A?+ i T. ix ) ir4s ‘ 
L!J Jr (s— s o y>\r— r„) 3/i 


CD 2) 

The second term on the right side of equation (D2) is zero 
because Ap=0 along the leading edge. Subtraction of 
equation (Dl) from equation (D2) then yields 


XX 


Apdrodsa 


at (s-8,y'>(r-r a W\ 


Inasmuch as the integral equations of reference 5 are 
formulated in terms of Au and are inverted by means of 
Abel’s integral equation, only solutions satisfying the Kutta 
condition will result therein. 





bw_ A/I f r Audrads, 

bx~ 8x[J J x (a — s«) 3/2 (r — r®) 8/ * 

For points in region II of a cancellation 
Thus, for the wing of figure 3 (a) , 


(D3) 

bw . 
^ a^=°- 


f* ds„ 

" r A udr 0 

Jo (s— s 0 ) 8/2 J 

riC.) 0 fa) 3/2 


(D4) 


This equation is the same as equation (11a) except that Au 
replaces A <p. The inversion by Abel’s integral equation for 
Au n in terms of A w r then gives (from equation (15)) 


Au, 


yr— r 2 (s) f r a (, > A utdr, 


r* j(«) 

JnW fa- 


ir— r 0 )Vr,(a)— r„ 


(24a) 





(s) Original position of wing. 

(b) Wing moved upstream distance ir. 

FlQiml 9.— Areas of Integration relating to equations (Dl) and (D2). 
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